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Abstract 

Consider a linear autonomous Hamiltonian system with a time periodic bound state solu- 
tion. In this paper we study the structural instability of this bound state relative to time 
almost periodic perturbations which are small, localized and Hamiltonian. This class of 
perturbations includes those whose time dependence is periodic, but encompasses a large 
class of those with finite (quasiperiodic) or infinitely many non-commensurate frequencies. 
Problems of the type considered arise in many areas of application including ionization 
physics and the propagation of light in optical fibers in the presence of defects. The mech- 
anism of instability is radiation damping due to resonant coupling of the bound state to 
the continuum modes by the time-dependent perturbation. This results in a transfer of 
energy from the discrete modes to the continuum. The rate of decay of solutions is slow 
and hence the decaying bound states can be viewed as metastable. These results generalize 
those of A. Soffer and M.I. Weinstein, who treated localized time-periodic perturbations of 
a particular form. In the present work, new analytical issues need to be addressed in view 
of (i) the presence of infinitely many frequencies which may resonate with the continuum 
as well as (ii) the possible accumulation of such resonances in the continuous spectrum. 
The theory is applied to a general class of Schrodinger operators. 
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1. Introduction 



1.1. Overview 

Consider a dynamical system of the form: 

id t (j> = H 0, (1.1) 

where H denotes a self-adjoint operator on a Hilbert space 7i. We further assume that H has 
only one eigenstate ipo G 7i with corresponding simple eigenvalue \$. Thus, 

&„(*) = e - lAo Vo (1.2) 



is a time-periodic bound state solution of the dynamical system (|1.1|) . We next introduce the 
perturbed dynamical system: 

idt<f> = ( H + eW(t) )<t>. (1.3) 

In this paper we prove that if the perturbation , eW(t), is a small, "generic" and almost periodic 
in time [], then solutions of the perturbed dynamical system (|1.3f) tend to zero as t — > ±oo. It 
follows that the state, &*(£), does not continue or deform to a time periodic or even time almost 
periodic state. Thus, 6*(t) is structurally unstable with respect to this class of perturbations. 
Our methods yield a detailed description of the transient (t large but finite) and long time 
(t — ► ±oo) behavior solutions to the initial value problem. Theorems |2.1[ - |2.3| contain precise 
statements of our main results. The following picture emerges concerning time evolution ( |1.3| ) 
for initial data given by the bound state, ipo, of the unperturbed problem. Let 

P(t) = |(^o,<K*))| 2 , (1.4) 

the modulus square of the projection of the solution at time t onto the state ipo 0- Then, 

(i) P(t) ~ 1-CV \t\ 2 , for |*| smallfj 

(ii) P{t) ~ exp{-2e 2 Tt) for t < O^r)- 1 ), T = 0{W 2 ), and 

(iii) P(t) ~ (t)~ a for |t| » (£ 2 r)- x , for some a > 0. 

1 See the appendix in section ^ as well as [0, || for definitions and results on almost periodic functions. 

2 (/, g) denotes the inner product of / and g. If ipo is normalized then P(t) has the quantum mechanical 
interpretation of the probability that the system at time t is in the state ipo. 

3 We do not discuss the short time behavior in this article; see fl^| . This small time behavior is related to 
the "watched pot" effect in quantum measurement theory p5[. 
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The time r = (e 2 r) _1 is called the lifetime of the state &*(£), which can be thought of 
as being metastable due to its slow decay. The mechanism for large time decay is resonant 
coupling of the bound state with continuous spectrum due to the time-dependent perturbation. 
Our analysis makes explicit the slow transfer of energy from the discrete to continuum modes, 
and the accompanying radiation of energy out of any compact set . 

Phenomena of the type considered here are of importance in many areas of theoretical physics 
and applications. Examples include: ionization physics |3|, f|, [TTJ and in the propagation of light 
in optical fibers in the presence of defects [O; see the discussion below. 



The results of this article generalize those of Soffer and Weinstein where the case: 

W{t) = cosfr/t) /3, (3 = (3* (1.5) 

was considered. The method used is a time dependent / dynamical systems approach introduced 
in [[H]], [^3] for the problem of perturbations of operators with embedded eigenvalues in their 



continuous spectra, [24], in the context of resonant radiation damping of nonlinear systems, as 



well as in p2 |; see also [|T^]. New analytical questions must be addressed in view of (i) the 



presence of infinitely many frequencies which may resonate with the continuum as well as (ii) 
the possible accumulation of such resonances in the continuous spectrum. This leads to a careful 
use of almost periodic properties of the perturbation (Theorems [2.1| and |2.2|) and hypothesis 
(H6) (Theorem |2.3|) , which is easily seen to hold when the perturbation, W(t), consists of a 
sum over finite number of frequencies, fij. 

A special case for which the hypotheses of our theorems are verified is the case of the 
Schrodinger operator Hq = —A + V(x). Here, V(x) is a real- valued function of x G H 3 
which decays sufficiently rapidly as \x\ — > oo. In this setting Soffer and Weinstein studied 
in detail the structural instability of Kit) by considering the perturbed dynamical system fli~3]), 
with W(x, t) = (3{x) cos(fit). Here, we consider a class of perturbations of the form W(x, t) = 
J2j Pj( x ) cos fijt, where the sum may be finite or infinite and where the frequencies fij need 
not be commensurate, e.g. W(x,t) = (3\{x) cost + ^{x) cosv^2t, where @i(x), i = 1, 2 is 
rapidly decaying as x — > oo. 

In addition to the problem of ionization by general time varying fields, we mention other 
motivations for considering the class of time dependent perturbations sketched above and defined 
in detail in section ||. 

(a) An area of application to which our analysis applies is the propagation of light through an 
optical fiber |13| . In the regime where backscattering can be neglected, the propagation of waves 
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down the length of the fiber is governed by a Schrodinger equation: 



id z (j> = ( -A ± + V(x±) )(j> + W(x ± ,z)<j>. (1.6) 

Here, <fi denotes the slowly varying envelope of the highly oscillatory electric field, a function 
of z, the direction of propagation along the fiber, and x± G IR 2 , the transverse variables. 
V(x±) denotes an unperturbed index of refraction profile and W(x±,z), the small fluctuations 
in refractive index along the fiber. These can arise due to defects introduced either accidentally 
or by design. The models considered allow for distributions of defects which are far more general 
than periodic. Our analysis addresses the simple situation of energy in a single transverse mode 
propagating and being radiated away due to coupling by defects to continuum modes. The 
bound state channel sees an effective damping. In particular the results of this paper have been 
applied to a study of structural instability of so-called breather modes of planar "soliton wave 



guides" [12]. The case of multiple transverse modes is of great interest [13|. Here one has the 
phenomena of coupling among discrete modes as well as the coupling of discrete to continuum 
/ radiation modes 0. There is extensive interesting work on this problem in the case where 
W(x±, z) is a stochastic process in z and radiation is neglected ||. 

(b) Nonlinear problems can be viewed as linear time- dependent potential problems where the 
time-dependent potential is given by the solution. A priori one knows little about the time 
dependence of the solution of a nonlinear problem. Nonlinearity is expected, in general, to 
excite infinitely many frequencies. Therefore results of a general nature for potentials with 
very general time dependence are of interest. This point of view is adopted by I.M. Sigal [jHJ, 
|2H|] , who considers the case where the nonlinear term defines a time-periodic perturbation, and 
then proceeds to study the resonance problem via time-independent Floquet analysis applied to 
the so-called Floquet Hamiltonian. The dilation analytic techniques used were first applied in 
the context of time-periodic Hamiltonians by Yajima [|^, |28f . Floquet type methods were 
also used in the time-periodic context by Vainberg ||25|| . The general class of perturbations we 



consider are not treatable by Floquet analysis and time-dependent analysis appears necessary. 
1.2. Outline of the method 

We now give a brief outline of our approach. For simplicity consider the initial value problem: 

id t (j)(t,x) = H (j)(t,x) + eW(t,x) (j)(t,x), (1.7) 
0|t=o = 0(0) (1.8) 
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where 

H = -A + V{x), W(t, x) = git) p(x), git) = £ 9j (1.9) 

j 

is a real- valued almost periodic function of t, and (3{x) is a real- valued and rapidly decaying 
function of x as \x\ — ► oo. The unperturbed problem (e = 0) can be trivally written as two 
decoupled equations governing the bound state amplitude, a(t), and dispersive components, 
4>d(t), of the solution. Specifically, let 

4>{t) = a(t)Mx) + <MM), (Vo,0d(*)) = 0. (1.10) 

Then, 

id t a(t) = X a(t), 
idt<f>d(t,x) = H (f) d (t,x), (1.11) 

with initial conditions: 

a(0) = (^o,0(O)), 
fa(P) = P C 0(O), (1.12) 

where 

defines the projection onto the continuous spectral part of H . 

For initial data a(0) = 1, 0^(0) = 0, we have a(t) = e~ tXot , <pd{t) = 0, corresponding to the 
bound state, b*(t). 

We now ask: 

(a) Under the small perturbation eW(t,x) does the bound state deform or continue to a nearby 
periodic or even almost periodic solution?, 

(b) How do solutions to the perturbed initial value problem behave as \t\ — > oo? 

For small perturbations sW(t,x) it is natural to use the decomposition (|1.10|) . Substitution 



of ( P-.10I) into ( p..3|) yields a weakly coupled system for a(t) and (pdif). This system is derived 
and analyzed in detail in sections 4-6. 

In order to illustrate the main idea, we introduce a simplified system having the same general 
character: 

id t a(t) = A a(t) + eg(t) (f3i/j , <fi d (t)) 
id t (f) d (t,x) = -A(f) d (t,x) + ea{t)g{t)P{x)%{x). (1.13) 
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Here, we have replaced Hq on its continuous spectral part by —A. 

If e(3 is small then A(t) = e lXot a(t) is slowly varying (dtA(t) = 0(ef3)). In particular, we 
have 

id t A(t) = ee lXat g(t) ( (3^oAd(t) ) 
id t (t> d {t,x) = -A(f) d (t,x) +A(t)e- iXot eg(t)f3(x)Mx). (1.14) 



Viewing A(t) as nearly constant, we see that the inhomogeneous source term in ( |1.14p has 



frequencies Ao + Pj] see (|1.9|). Therefore, if Ao + Pj > 0, for some j then Ao + fij lies in the 
continuous spectrum of —A (Hq) and therefore <frd satisfies a resonantly forced wave equation. 
A careful expansion and analysis to second order in the perturbation eW(t) (see the proof of 
Proposition 4.1) reveals the system for A(t) and 4> d (t) can be rewritten in the following form, in 
which the effect of this resonance is made explicit: 

d t A(t) = (-e 2 T + p(t) ) A(t) + E(t;A(t),<p d (t))- (1-15) 
idMt.x) = H 4> d {t,x) + P c F(t,x;A(t),(f) d (t)). (1.16) 

The terms E(t) and F(t, x) formally tend to zero if A(t) tends to zero and if the "local energy" 
of 4> d {t) tends to zero as t — > oo. The strategy of sections 5 and 6 is to derive coupled estimates 
for A(t) and a measure of the local energy of <p d from which one can conclude, for eW(t) small, 
that solutions to ( |1.15| - |1.16|) decay in an appropriate sense. The key to the decay of solutions is 
the constant T, given by 

r = - E \mf (PcWo, ^o-Ao-^Pe^o); (i.i7) 

{j : A o +Mj>0} 

see also hypothesis (H5) of section 0. The quantity T is a generalization of the well known Fermi 
golden rule arising in the theory of radiative transitions in quantum mechanics ||, [|, [l(| . For the 



example at hand, Q1.9D , the sum in ( |1.17| ) is over all j for which /ij + Ao is strictly positive, i.e. lies 
in the continuous spectrum of Hq. Thinking of Hq as having a spectral decomposition in terms 
of eigenfunctions and generalized eigenfunctions, let e(A) denote a generalized eigenfunction 
associated with the energy A. Then each term in the sum (|1 . 1 7|) is of the form: 

|(e(A + ^),Wo)| 2 , (1-18) 

Thus clearly T > 0, generically. 

Neglecting for the moment the oscillatory function p(t) in ( |1.15| ), we see that coupling of the 
bound state by the time dependent perturbation to the continuum- - radiation modes, at the 
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frequencies pj + Ao > 0, leads to decay of the bound state. The leading order of equation ( |1.15 
|1.16|) is normal form in which this internal damping effect is made explicit; energy is transferred 
from the discrete to the continuous spectral components of the solution while the total energy 
remains independent of time: 

II <t>(t) 111 = Ht)\ 2 + || <f> d (t) \\l 

= |a(0)| 2 + || 0,(0) \\l (1.19) 



1.3. Energy flow; contrast with the analysis of [22 



The goal is to show that energy flows out of the bound state channel into dispersive spectral 
components. The normal form above is the system in which this energy flow is made explicit. 
Once the normal form (|1.15| - |lTT6|) has been derived, it is natural to seek coupled estimates for 
A{t) and <fid{t) from which their decay can be deduced. This is implemented in section |(| A 
natural first step is to introduce the auxiliary function: 

A{t) = efi ds A(t), (1.20) 

for then A(t) satisfies simplified equation of the form: 

d t A(t) = -s 2 TA(t) + E(t;A(t),Mt)) (1-21) 

If 3ft Jq p(s) ds is uniformly bounded then, modulo time-decay estimates on E(t; A, <pd) and 
F(t;A, (f)A, the decay of A(t) and therefore of A(t) follows. For the class of perturbations 
considered in |22j p(t) is a periodic function, having only a finite number of commensurate 
frequencies, none of them zero. Therefore, in this case 3ft Jq p(s) ds is uniformly bounded. 
However, in the present case p(t) is almost periodic with mean M(3ftp) = (see section |9[); p(t) 



is displayed in (^4.12|) . 3ftp(t) has, in general, infinitely many frequencies, pk — Pj, k ^ j which 
may accumulate at zero. Most delicate is the case where, along some subsequence, pk — Pj — > 0. 
It is well known that the integral of an almost periodic function of mean zero is not necessarily 
bounded 0, so we are in need of a strategy for estimating the effects of 3ft J * p(s) ds. We 
address the estimation of 3? J * p(s) ds in two different ways corresponding to Theorem 
(section |5.1| ) and Theorems p7T] - pT2] (section [5 . 2| ) . In section |5]1| 3? J * p(s) ds is estimated under 
the hypothesis (H6) which requires that the rate of accumulation of a subset of frequencies 
{pj}j £ i is balanced by the decay of the Fourier coefficients gj as j — > oo, j e /. This leads to 
a bound on 3ft Jq p(s) ds (Proposition |5.2| ). In section f>.2\ the estimates are based on a more 
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refined analysis; the almost periodic function pit), is decomposed into a part with bounded 
integral and a part which has mean zero. The latter is controlled using results on the rate at 
which an almost periodic function approaches its mean. 



1.4. Fermi golden rule and obstructions to Poincare continuation 

In the theory or ordinary differential equations it is a standard procedure, given a periodic 
solution of an unperturbed problem, to seek a periodic or almost periodic solution of a slightly 
perturbed dynamical system. We now investigate this procedure in the context of Ql.Tf) and its 
solution b*(t) for e — 0. Seek a solution of the form: 

<f>(t) = K(t) + (f>i{t) + 0{e 2 (3 2 ). (1.22) 

Here, fa = 0(e/3)£j Substitution of (|T22D into ([LTD yields the equation: 

idtfc = H fa + e/3 g(t) b*(t). (1.23) 

This equation has a solution in the class of almost periodic solutions of t with values in the 
Hilbert space TC only if (3 g(t) b*(t) is "orthogonal" to the null space of id t — H . 

We now derive this condition. Let e(£) be a solution of H e(() = C e (0- Then, taking the 
scalar product of (|1.23| ) with e - ^* e(() and applying the operator lim^oo T _1 J T • dt to the 



resulting equation gives: 

= lim T- 1 / e* 1 e- lXot g(t) dt ( e(C),^o ) • (1-24) 

^Too Jo 

Substitution of the expansion for g(t) yields: 

gjS(C,X + ^) (e(C),/3Vo) = 0, (1.25) 

j'GZ 

where 5(a,b) = if a ^ b and 5(a, a) = 1. If C > which lies in the spectrum of H , satisfies 
C = Ao + yUfc for some k E 71 (which will be the case in our example if Ao + yUfc > 0), then we 
have that: 

(e(Ao + /i fc ),/fyo) = (1.26) 



is a necessary condition for the existence of a family of solutions of (E7) which converges to 

as the perturbation tends to zero. We immediately recognize the inner product in ( |1.26| ) as 

the projection of j3ipo onto the generalized eigenmode at the resonant frequency Aq + yUfc, which 



arises in ( |1.17| ); see also ( |1.18| ). Therefore the obstruction to continuation of K(t) to a nearby 



almost periodic state of the system can be identified with the damping mechanism. 



J This argument is heuristic so we do not specify the norm with which the size of f3 is measured. 
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1.5. Outline 



The paper is structured as follows. In section |2| we give a general formulation of the problem. 
The hypotheses on H , the unperturbed Hamiltonian and W(t), the perturbation are introduced 
and discussed. There are two types of theorems: Theorems |2.1| & |2.2j and Theorem |2.3| . Although 
the conclusions of these are quite similar, as discussed above, they differ in a key hypothesis on 
the perturbation W(t), which is relevant in the case where W(t) has infinitely many frequencies 
which may resonate with the continuous spectrum. In section || we apply the results of section 
| to the case of Schrodinger operators H = —A + V(x) defined on L 2 (H 3 ). To check 
the key local energy decay hypotheses we use results of Jensen and Kato || on expansions 
of the resolvent of Ho near zero energy, the edge of the continuous spectrum. In section [| 
the dynamical system (|1.3| ) is reformulated as a system governing the interaction of the bound 
state, and dispersive part of the solution. This section contains an important computation, in 
which the key resonance is made explicit and a perturbed "normal form" for the bound state 
evolution is derived (Proposition 4.1). Sections |5| and ^| contain estimates for the bound state 
and dispersive parts of the solution for intermediate and large time scales. In section [7] we 
discuss extensions of our Theorems [2.1| - [2.3| to a more general class of perturbations. We shall 
frequently make use of some singular operators which are rigorously defined in section |8], an 
appendix, and of elements of the theory of almost periodic functions || M, which are assembled 
in section |S], the second appendix. 

Notations and terminology: 

Throughout this paper we will use the following notations: 

IN ={1, 2, 3,...}; 

IN ={0, 1, 2, 3,...}; 

TL ={...,-3, -2, -1, 0, 1, 2, 3,...}; 

for z a complex number, ^kz and ^sz denote, respectively, its real and imaginary parts; 
a generic constant will be denoted by C, D, etc; 
(x) = (l + |x| 2 )^; 

C(A, B) = the space of bounded linear operators from A to B; C(A, A) = C(A). 

Functions of self-adjoint operators are defined via the spectral theorem; see for example [IT? 



The operators containing boundary value of resolvents or singular distributions applied to self- 
adjoint operators are defined in section ||. 

Acknowledgements: This research was supported in part by National Science Foundation 
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2. General formulation and main results. 

Consider the general system 

id t (j)(t) = (H + W(t))<f>(t), 
<f>\t=o = 0(0). (2.1) 

Here, <fi(t) denotes a function of time, t, with values in a complex Hilbert space 7i. 
Hypotheses on H : 

(HI) H is self-adjoint on H. and both H and W(t), t G IR 1 , are densely defined on a subspace 
V of U. 

The norm on TL is denoted by || • ||, and the inner product of /, g G H, by (/, g). 
(H2) The spectrum of H is assumed to consist of an absolutely continuous part, cr cont (H ), 
with associated spectral projection P c and a single isolated eigenvalue Ao with corresponding 
normalized eigenstate, ipo, i-e- 

H ij = A ^o, HVoll = 1- (2-2) 

The manner in which we shall measure the decay of solutions is typically in a local decay 
sense, e.g. for the scalar Schrodinger equation governing a function defined on H n we measure 
local decay using the norms: / i— > || {x)~ s f\\i J 2 ) where s > 0. So that our theory applies to a 
class of general systems (involving, for example, vector equations with matrix operators), we 
assume the existence of self-adjoint "weights", W- and w + such that 

(i) w + is defined on a dense subspace of 7i and on which w + > cl, c > 0. 

(ii) w_ G C(TC) such that Range(w _) C Domain{w + ). 

(iii) w + W- P c = P c onH and P c = P c W- w + on the domain of w + . 

In the scalar case w + and w_ correspond to multiplication by (x) s and (x)~ s , respectively, 
see section |3|. 

The following hypothesis ensures that the unperturbed dynamics satisfies sufficiently strong 
dispersive time-decay estimates. Let {/i J } j . gZ denote the set of Fourier exponents associated 
with the perturbation W (see hypothesis (H4) below). 
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(H3) Local decay estimates on e lHot : Let r\ > 1. There exist w + and as above, and a 
constant C such that for all / e H satisfying u> + / EH we have: 

(a) || W _ e -^*P c /|| < C^-^IK/H, forte H; (2.3) 
(b) ||u;_ e - <ff ° t (Ho-Ao-/x i -iO)- 1 Pe/|| < C (t)- ri |K/||, fort>0 (2.4) 



and for all j E TL. For t < estimate (|2.4j) is assumed to hold with — iO replaced by -MO. See 
section |8| for the definition of the singular operator in (|2.4|) 



Remark 2.1. There is a good deal of literature on local energy decay estimates of the form 
form for e - iHot P c in the case H = —A + V(x) on L 2 (H n ). These results require sufficient 
regularity and decay of the potential V(x). We refer the reader to jj^j, J^/ and fT^/; see also jj^j , 



Remark 2.2. Estimates of the type (H3b) are obtained in t 22, |23|, Appendix A]. A key point 
here is that we require that one can choose the constant, C, in Qp.4\) to hold for all fij. It appears 
difficult to deduce this uniformity of the constant by the general arguments used in and 
However, in section where we apply our results to a class of Schrddinger operators, we 
can verify (H3b) using known results on the spectral measure. 

(H4) Hypotheses on the perturbation W(t): 

We consider time-dependent symmetric perturbations of the form 

W{t) = -fa + cos(^t) /3j, with /3* = fa and £ \\Pj\\c(H) < oo. (2.5) 

In many applications, f3j are spatially localized scalar or matrix function. Note that formula 
(|2.5|) can be rewritten in the form: 

W{t) = \ ^exp(-^t)/3 i , (2.6) 

where, /i = and for j < 0, fij — —fi-j, Pj = P-y Thus, W(t) is an almost periodic function 
with values in the Banach space £(7~C), with the Fourier exponents {A*jLj gZ and corresponding 
Fourier coefficients {/3j}- EZ ; see, for example, ||. 

To measure the size of the perturbation W, we introduce the norm 

infill = -EIK&IUpo + \\w-,n + ), (2.7) 

which is assumed to be finite. Here TC + , respectively 7Y_, denote the closure of the domain of 
w + , respectively the range of P c , with norm / — > respectively / — > 
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Remark 2.3. A special case which arises in various models, is: 



W{t) = g(t)P, (2.8) 

where 

g(t) = J2 9j cos fijt, (2.9) 

j 

\\w + j3\\c(H) + 1 1/3 1 1 < °o and the sequence {gj} is absolutely summable. 

Remark 2.4. Our results are valid in the more general case 

Wit) = -p +J2 cos^t + 5,) fy, 
z jew 



where Pj are self-adjoint such that expression ( j2.7p is finite. This follows because the proofs use 



only the self-adjointness ofW and the expansion: 

W ( f ) = o J2 expi-tfijt)^, 

where = e^Wifc and ^ = -fij, ^ = 0. 

We will impose a resonance condition which says that {A + H cr cont (-£fo) is nonempty 

and that there is nontrivial coupling; see section 1.4. Let us first denote by I res the following 
set: 

I res = {j e 7L : X + Hj e a cont (H )}. (2.10) 

(H5) Resonance condition - Fermi golden rule I res is nonempty and furthermore, there exists 
#o > 0, independent of W, such that 

r = \ E (*Wo, 5(if - A - //^PcAVo) >e |||W||| 2 >0 (2.11) 

jdzlres 



Remark 2.5. For the exact definition of the Dirac type operator in ( \2. 1 1\ ) see section [|. That 
T is finite is a consequence of the estimate ( |8.8| ) and 



T <^J2\\™+Pj\\ 2 <—\\\W\\\ 2 (2.12) 



7T 7T 

3 



We now state our main results: 
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Theorem 2.1. Let us fix Hq and W(t) satisfying hypotheses (H1)-(H5). Consider the initial 
value problem: 

id t( j)(t) = {H + eW(t))<f>(t), 

0U=o = 0(0), (2.13) 

with w + (f>(0) G 7i. Then, there exists an Eq > (depending on C, r\, and 9q) such that whenever 



\e\ < e , the solution, <p(t), of ( \2.1 3j ) satisfies the local decay estimate: 

Ik- <f>(t)\\ < C{t}- ri \\w + 0(O)||, t G H. 



(2.14) 



Under the same hypotheses as Theorem [2.1| , we obtain more detailed information on the 
behavior of <f)(t): 



Theorem 2.2. Assume the hypotheses of Theorem 2Jj. For any < 7 < T there exist the 



constants C and D (depending on C, r±, 9 and ~f) such that any solution of ( |2.13| ), for \e\ < e 
and w + 4>(0) G 7i, satisfies: 



(f){x,t) = a(t)ipo + (f) d (t), (V>o , (j>d{t)) 



0. 



a(t) 
P(t) 



a(0) e 
P(0) e 



-£ 2 (r-7)l*L^W 



+ Ra{t) 



-2 £ 2 (r- 7 )|*| 



+ Kit) 



-iH t 



P C 0(O) + 4>(t). 



(2.15) 



where T is given by l\2.11\ ) and u)(t) is a real-valued phase given by 
co(t) = X Q t-e (Vo, J W(s)ds ^0) 



+ -A E (Prto, p.v.(if - Ao - ^r'PcP^ 



P(t) is defined in ( \1.^ ) and for any fixed T > we have 



\R a (t)\ < C \e\ \\\W\\\, \t\ < 



(2.16) 







T 




< 




e 2 T 






T 


t\ 


< 




e 2 T 



(2.17) 
(2.18) 



Moreover, 



\R a (t)\ = wr 1 ), = o((t)-^), \t\ -> 00. 
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Finally, = 0i + 02 is given in (|4.fip, with \\w-(f>(t)\\ 
(H3) same holds \\w-<p d (t)\\ = 0((t)- ri ) as \t\ oo. 



0({t) ri ) as \t\ — > oo. Therefore, by 



Remark 2.6. Suppose the initial data is given by the bound state of the unperturbed problem, 
i.e. 0(x, 0) = ipo(x), a(0) = 1, 0d(O) = 0. Then, from the expansion of the solution we 
have that for < t < e _2 r _1 that P(t) (see is of order e~ 2e2(T ~ 7 )*, with an error of 

order e. Hence it is natural to view the state ip e~ lXot as a metastable state, with lifetime 



T 



\W\ 



Although 7 > is arbitrary we have not inferred that the 



actual lifetime is r = e^T^ 1 under hypothesis (H1)-(H5). The reason is that the constants C 
and D in the estimates ( |2. 1 7\ ) and ( ^.18j) blow up as 7 \ 0. In order to remedy this we need an 
additional hypothesis: 



(H6) Control of small denominators : There exists £ > 0, independent of W, such that 

1 



E 



fij — fi k 



<Z\\W 



(2.19) 



Remark 2.7. By ( 18.51) we have 



£ \{M ,5{H Q -\ - N )^)\<C t:- 1 \\\W\ 

jei res , feez, j+k 



is finite (see also Remark |2. 5p . Thus, (H6) is important only if: 

mi{\nj -fi k \: j,k G j ^ k and A + fij E cr cont (H )} = 0, 



(2.20) 



(2.21) 



i.e. the Fourier exponents {fij} are such that Ao + fij accumulate in a c . In particular, if the 
perturbation W(t) consists of a trigonometric polynomial: 



N 



W(t) = Y cos fj,jt Pj, 



(2.22) 



then (H6) is trivially satisfied. 



Remark 2.8. Hypothesis (H6) can be imposed by balancing the clustering of the frequencies 
A + fij in the continuous spectrum of H with rapid decay of (Pkipo , S(H — A — fij) Pjipo) as 
j, k — > 00. Let (3j(x) = Qj (3(x). Then, W(t,x) = J2j (7j cos(//_ 7 t) f3(x). Using Remark |2T5| we 
find that the left hand side of (\2.19j ) is bounded by Y.j.k£Z-.j±k \9j9k\ \fij — A^l -1 lll^lll 2 - The 
constant £ in ( |2. JS| ) is finite if, for example, fij = 2\X \ + Qj = |j|~ 2 ~ r , t > 0. 
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In case (H6) is satisfied we have the following improvement of Theorem |2.2| : 



Theorem 2.3. Assume the hypotheses (H1)-(H6) hold. Then there exist Eq and the constants 



C, D (depending on C, r±, 8 and such that any solution of ( |2J3[ ), for \e\ < eq and w+(f)(0) G 
TC, satisfies 

<f){x,t) = a(t)^ + (f) d (t), (z/> , <Pd(t)) = 0, 

a(t) = a(0) e™^ 1 ' 1 e luj{t) + R a (t) 

P(t) = P(0) e~ 2£2m + R' a (t) 

Mt) = e~ lHot P C 0(O) + 4>{t). (2.23) 



Here, u(t) is given by ( \2.1Q ) and R a (t), R' a (t), W-(f>d(t) satisfy the estimates of Theorem |3Q . 



3. An application: the Schrodinger equation 

In this section we verify hypotheses (H1)-(H4) in the particular case of the Schrodinger equation 
on the three dimensional space with a time almost periodic and spatially localized perturbing 
potential: 

id t <p = (-A + V(x))(f) + £W(x,t)(f>, (3.1) 
with : IR 3 x IR — > C, (x, t) — > 0(x, t) and 

W(x,t) = ^(3 (x) + cos^t)^-^) 

where fa G IR, j G HST , and (3j : IR 3 — > IR, j G IN, are localized functions. Models of the sort 
considered in this example occur in the study of ionization of an atom by a time- varying electric 
field; see [0, §. 



We take TC = L 2 (IR 3 ), and H = —A + V(x), where V(x) is real- valued with moderately 
short range. More precisely, we suppose that there exists a > 4 and a constant D such that 

\V(x)\ <D(l + \x\)- a . (3.2) 

Thus, Hq is self-adjoint and densely defined in I? . In what follows we assume that Ho has exactly 
one eigenvalue which is strictly negative and that the remainder of the spectrum is absolutely 
continuous and equal to the positive half-line. Our results can be extended to operators with 
strictly negative, multiple eigenvalues 0. 
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We first discuss the local decay hypothesis (H3). As weights used to measure local energy 
decay we take w± = (x) ±s where s > 7/2 and fix ri = 3/2. Our aim is to obtain the estimates: 

(H3a) Ww^e-^PJW < C (t)-^\\w + f\\, (3.3) 
(H3b) IKe-^ifo-Ao-^-iO^Pc/ll < C (t)^ 2 \\w+f\\, (3.4) 

for all fij G /Z, with C independent of j. 

We shall assume that the frequencies {A + yUj} do not accumulate at zero, the edge of the 
continuous spectrum of H : 

m* = min{ | Aq + Mi I : 3 ' e 2Z } > 0. (3.5) 



To prove (|3.3| ) and (|3.4j ) we use the spectral representation for the operators e tHot P c and 



e -iffot ( Hq _ A - fij - iO) _1 P c , namely: 

e -*v tp c = r e ^t E >( X )d\ (3.6) 







poo 

iHot {H - A - m - iO) _1 P c = / e- iA '(A-A -^ -iO)- 1 E'(A)c/A (3.7) 



where E'(X) = n^ 1 ^(H — A — iO)^ 1 is the spectral density induced by H , ||. 

The technique of getting (H3a) from ( |3.6|) is presented in || Section 10] and it can be 
summarized in the following way. We decompose the integral in (|3.6| ) in two parts, corresponding 
to low energies (A near zero) and high energies (A away from zero) by writing 

E > = X E' + (1- X )E' 

roo roc 

w _ e ~iH tp cW _ = / e -*^x) w _E'(X)w-dX + / e~ iXt (l - x(X))w-E'{\)w-d\ 
Jo Jo 
= S 1 + S 2 . (3.8) 

Here xW * s a smoothed characteristic function of a neighborhood of origin, chosen so that 

X(A) = 1, |A| < ^m* 

3 

X(A) = 0, |A| > -m*. 

To estimate the two integrals in ( |3.8| ) we make use of the detailed results of on the family 
of operators { -E"(A) }. First, by Theorem 8.1 and Corollary 8.2 of 0, W-d^E'(\)w- is bounded 
on L 2 and satisfies 

|| W-d\E'{\)w- || B(L 2) = 0(A- (fe+1)/2 ) as A -> oo, (3.9) 
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for k G {0,1,2,3}. Integration by parts twice in the second integral in Q3.8|) and use of the 
estimate ( |3.9|) with k = 2 yields the estimate: 

II S 2 \\b(l 2 ) — °(t~ 2 ) as t — > oo. (3.10) 

Next, by Theorem 6.3 of || we have the low energy asymptotic expansion: 

W-E'(X)w- = -A- 1/2 5_i + \ l/2 B 1 + o(A 1/2 ) as A -> 0, (3.11) 

where f?_i, Bi are bounded linear operators on L 2 . Use of this expansion in the first integral 
of (|3.8|) yields the expansion in B(L 2 ): 

Si = (7Ti)- 1 /2 r l/2 jB _ i _ ( 47r ,)-l/2 r 3/2 5i + ( f -3/2) ag t _> ^ ( 3 . 12 ) 

Thus (H3a) is satisfied provided that -B_i is the null operator or equivalently H^ip = has 
no solution with the property W-ip G L 2 (IR 3 ). The last condition holds for generic potentials 
V(x) and when it is violated one says that H has zero energy resonance; see for details. 

In the same way one can prove (H3b) from the spectral representation fl3.7D provided that 
the integral is non-singular, i.e. Ao + fij < 0. In the case Ao + yUj > m* > we first decompose 
the singular integral in two parts, one away from singularity point, Ao + fij, and the other in a 
neighborhood of it by using the smoothed characteristic function 

XjW = X(A-A (3.13) 
which is supported in a neighborhood of Ao + fJ>j, which does not include A = 0: 

POO 

e -*Hot( HQ _ Xo _ tlj _ iQ yip e = / e-^A-Ao-^-^l-x^A))^)^ 

J 

roc 

+ / e- iXt {X- \ - fij-iO^xjiXjE'iXjdX, (3.14) 

J 

The non-singular integral may be treated as above while the singular one defines the singular 
operator: 

Tj = e- lHot (H - Ao - m - zO)- 1 Xj (^o)Pc, 
via the spectral theorem. Here, Tj = ]xm v \oTj , where 

TV = e- tHot (H - Ao - m - i^x^V,. 

To estimate its L 2 operator norm we use the integral representation 

w_T]w_ = - / e^^^^w^e-^xAH^cW-ds. (3.15) 
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But this reduces to the evaluation of 

roo 

w„e~ tHoS Xj(H )P c w_ = / e- lXs Xj (\)w_E'(\)w_d\, s>t, (3.16) 

J 

where we used again the spectral representation theorem. Integration by parts three times in 
Q3.16| ) and use of the estimate ( |3.9D with k = 3 implies 

\\w-e' iHoS Xj(H )P c W-\\ B{L 2 ) = o(s~ 3 ) as t -> oo. 

Replacing this in (|3.15|) , integrating and passing to the limit as r\ \ we obtain an o(t~ 2 ) 



estimate for Tj which is even better than we need to satisfy (H3b). 

Moving now towards hypothesis (H4), we may choose the time-dependent perturbation to 
be of the form: 

W{x,t) = \(3 + ^cos^^(x), (3.17) 

with Pj rapidly decaying in x, e.g. (x) 2s \\{3j(x) \\ < Cj for all x e M 3 , j E 1N , where SjgN Q < 
oo. Thus, (H4) is satisfied as well. 



Therefore, our main results Theorems 2.1-2.2 on the structural instability of the unperturbed 



bound state, and large time behavior for systems of the form (3T) apply provided (H5), the 
Fermi Golden Rule resonance condition, holds. For results concerning more general perturba- 
tions than the ones in ( |3.1| ) see section 0. 

4. Decomposition and derivation of the dispersive normal form 

The results of this section rely on hypothesis (HI) through (H4) only, so they may and will be 
used in proving Theorems [2.1| - [2.3| . 

As in [21]], |2^] and |23[, we begin by deriving a decomposition of the solution, <p(t), which 
will facilitate the study of its large time behavior. Let 

<f>(t) = a(t)i> Q + (f> d (t), (4.1) 

with the orthogonality condition 

(V>o, &(*)) = fora11 *■ ( 4 - 2 ) 

Note therefore that (fid = P c (fid- 

We proceed by first inserting ( f4.1| ) into fl2.13| ), which yields the equation: 



id t a(t)ip + id t <t>d(t) = A a(t)^o + H <j> d {t) 

+ ea(t)W(t)ifi + eW{t)<p d {t) (4.3) 
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Taking the inner product of ( |4.3|) with ipo we get the following equation for a(t): 

id t a = X a(t) + e(ip ,W(t)i> )a(t) + e (V> , W{t)<j) d ) , (4.4) 
a(0) = (Vo,0(O)) 

In deriving ( |4.4|) we have used that i/jq is normalized and the relation 

(V>o,9^d) = 0, (4.5) 



a consequence of ( [4.2| ). 

Applying P c to (|4.3|), we obtain an equation for <p d : 



id t Mt) = HoMt) + eP c W{t)<f> d {t) + ea{t)P c W{t)ip , (4.6) 
0,(0) = P C 0(O) 

Since we are after a slow resonant decay phenomenon, it will prove advantageous to extract the 
fast oscillatory behavior of a(t). We therefore define: 

A(t) = e iXot a(t). (4.7) 

Then, ( |4.4[) reads 

d t A = -ieA(^ ,W(t)^ ) - tee^ (V>„, W(t) </>&)) ■ ( 4 -8) 
Solving (|4.6[) by Duhamel's formula we have 

Mt) = e~ mot M0) ~ ie e' tHo{t - s) P c W(s)a(s)^ ds 

Jo 

- ie f* e- iHo{t - s) P c W (s)<j> d (s) ds 
Jo 

= Mt) + Mt) + Mt)- (4.9) 



By standard methods, the system (gj)-(gj) for A(t) and Mt) = M) ~ e ~ iX ° l A (t) has a 
global solution in t with 

A e C\M), \\Mt)\\ e C°(1R), \\w^Mt)\\ e C°(IR). 

Our analysis of the |i| — > oo behaviour is based on a study of this system. 
By inserting Q4.9|) into Q4.8|) we get 

2 

$A(t) = -ieA(t){ik,W{t)iM - iee tXot Y,{^W{t)M- (4.10) 

3=0 
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We next give a detailed expansion of the sum in ( |4.10| ). It is in the j = 1 term that the key 
resonance is found. This makes it possible to find a normal form for ( |4.10| ) in which internal 
damping in the system is made explicit. This damping reflects the transfer of energy from the 
discrete to continuum modes of the system and the associated radiative decay of solutions. 



Proposition 4.1. Fort > 0, 

d t A{t) = ( -e 2 T + p(t) ) A(t) + E(t) 

where T is defined in (\2.11\ ), 

pit) = -ie{i/> , W(t)iM + 

+ -/ E (fl^o, P.V.(F - Ao - p^P^o 



(4.11) 



J2 e^-^!,. 



Vo, (H - A - fij - 20)~ 1 P c /3^o) (4.12) 



and 

m 



_L £ 2 e iX t 

4 



E ^ kt ( Mo, 



-e 2 A{0)e iXot £ ei 

j.A-GZ 

-iee iXot (^o, W(t)Mt)) 
-tee iXot (</>o, WW 2 (t)) 



Here, (fio and 4>2 are given in (\4.9j) . 



fc ^ , e~ iHot (Hq - Ao - fij - «0)- 1 P c /?^o 
■4W->)(H _ a - ^ - z0)- 1 P c e- i(Ao+ ^ )s 9 s A(s)/5 i ^ o j 

(4.13) 



Although the proposition is stated for t > 0, an analogous proposition with — e 2 T replaced by 
e 2 T holds for t < 0. The modification required to treat t < is indicated in the proof. 



Remark 4.1. (1) The point of (4. 11) is that the source of damping, T > 0, which arises due 
to the coupling of the discrete bound state to the continuum modes by the almost periodic 
perturbation, is made explicit. Note that dtp(t) is of order e 2 \\\W\\\ 2 as the first two terms of 
p(t) are pure imaginary inducing only a phase shift in the solution, A(t), while the last one is 
of the same order as the damping, and may compete with it. A key point of our analysis is to 
assess the contribution of this last term in ( 4.12^) . 

(2) The leading order part of equation ( |4. 11\ ) is the analogue of the dispersive normal form 
derived in pi/ for a class of nonlinear dispersive wave equations. 
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Proof of Proposition 4.1 

Using the expression for W(t) in (|2.6| ), which is a uniform convergent series with respect to 
t eE, and the definition A(t) = e iXot a(t), we get from ( p~9D 



0i (*) 



/ e- iHo{t - s) e- iXoS A(s)P c ]T e~^ s /3^o ds 



2 Jo 



is 



-iHo(t-s) -i(Ao+Ay)s 



(4.14) 



We would like to integrate by parts each of the integrals in the above sum. We cannot proceed 
directly since the resolvents of H in A + fij, j G 7L, would appear and hypothesis (H5) implies 
that some of the A j G 2L, are in the spectrum of H . Instead we regularize 4>i by defining: 



«(«)= - ? Ey < 



-iifo(*-*) e -*(^o+A 1 j+*^)s 



(4.15) 



for rj positive and arbitrary and t > 0. Note that 0i(t) = lim^\o0i(O uniformly with respect 
to t on compact intervals. 



Now, integration by parts for each integral in expression ( 4.15 ) and letting r\ tend to zero 
from above gives the following expansion of (ipo, W(t)(j)i(i)) : 



(Vo, W(t)Mt)) = W{t)ik, -- e~ iXot e- l ^A(t){H - A - ^ - zO)- 1 P c /3^o 



./'OZ 



+ ( ^Wo, ^(0) £ e" iffo *(^o - A - - z0)- 1 P c /3^o 
2 jez 



(4.16) 



+ (w(t)fo, |E j\- lH ^- s \H - A - //,• - 20)- 1 P c e- J ( A « + ^) s 5 s A( S )/?^o^j . 

The definition of the singular operators in the above computation is given in section [|[ The 
choice of regularization, +ir), in ( 4.15Q ensures that the latter two terms in the expansion of (f)i, 
( [4.16|) , decay dispersively as t — > +oo; see hypothesis (H3) and section §. For £ < 0, we replace 
with —27] in ( 4.15| ). 

To further expand the first series in ( |4.16|) we use the identities (|8.5|) . The proof of Proposition 
4~T| is now completed by substitution of (|8.5|) in the expansion (|4.16|) for X and of the result 



into the second term of the sum in ( 4.10 ). j 

In the next sections we estimate the remainder terms in (4.9) and ( |4.11| ). 
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5. Estimates on the bound state amplitude 



Our strategy is as follows. Equations ( f4.9| ) and ( |4.11| ) comprise a dynamical system governing 
4>d{t) and a(t) = A(t)e~ lXot , the solution of which is equivalent to the original equation (1.1). 
In this and in the following section we derive a coupled system of estimates for A(t) and <fid{t)- 
This section is focused on obtaining estimates for the bound state amplitude A(t) in terms of 
4>d{t), while the following section is focused on obtaining dispersive estimates for 4>d{t) in terms 
of A(t). We treat only the case t > since the modifications for the case t < are obvious. 
The coupled system of estimates shows that A(t) decays in time, provided 4>d(f) is dispersively 
decaying and vice- versa. We exploit the assumed smallness of the perturbation eW to "close" 
the resulting inequalities, and prove the decay of both A(t) and (j)d{t). 

The main difference from the strategy employed in |22] for the estimation of the bound state 
amplitude is related to the presence of infinitely many frequencies in the perturbation W(t). In 
particular one can have an accumulation of resonances in the continuous spectrum of H . We 
have two strategies for obtaining estimates for A(t) which correspond to the use of hypotheses 



(H1)-(H5) (Theorems 2A and |2.2| ) or hypotheses (H1)-(H6) (Theorem |2.3|) . These strategies 
revolve around estimation of 3? J * p(s) ds, where p is given by Q4.12|) . (H6), which controls 
certain "small divisors" which arise from the clustering of frequencies, ensures that 

& / p(s) ds < C £ 2 |||W||| 2 . (5.1) 
Jo 

This, in turn, implies that the contribution of p(t) in the size of A(t) is of order e 2 |||W||| 2 . 
Without hypothesis (H6) we carefully decompose p(t) as 

p(t) = e 2 a(t) + V (t), 

where a it) is a real almost periodic function with mean, M(o~), zero and 3? f r)(s) ds < 
Ce 2 \ \W 1 1 1 2 . As in the previous case, the contribution of the r)(t) in the size of A(t) is of order 
£ 2 | I IV^I 1 1 2 - On the other hand, a(t) competes with the damping term e 2 T in equation (|4.11|) , 
but being oscillatory (i.e of mean zero) and of the same size as the damping it allows the latter 
to eventually dominate. 

As the above discussion suggests it is simplest to start by assuming (H6) to get sharper 
estimates on A(t) (Theorem |2.3|) and then to relax this assumption (Theorem |2.2| ). We begin 
with a simple Lemma which we shall use in a number of places in this and in the next section. 

Lemma 5.1. Let a > 1. 

"' (t - s)~ a (s)~P ds < C a3 (5.2) 
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Proof: The bound is obtained by viewing the integral as decomposed into a part over [0, t/2] and 
the part over [t/2, t]. We estimate the integral over [0, t/2] by bounding (t — s)~ a by its value 
at t/2 and explicitly computing the remaining integral. The integral over [t/2,t] is computed 
by bounding (s)~@ by its value at t/2 and again computing explicitly the remaining integral. 
Putting the two estimates together yields the lemma. 

We now turn to the estimate for A{t) in terms of the dispersive norm of <pd{t) and local decay 
estimates for e~ lHot P C (H ) . 



5.1. Estimates for A{t) under the hypotheses of Theorem [2.3 



Proposition 5.1. Suppose (H1)-(H6) hold. Then A(t), the solution of $4.1 can be ex- 
panded as: 



A(t) 
R A (t) 



e IoP^ ds (e- e2rt A(0) + R A (t) 



(5.3) 
(5.4) 



where E(t) is given in ( |4.J5|) and ( |57^) . For any a > 1, there exists a 5 > such that R A (t) 
satisfies the estimates for T > 2(e 2 T)^ a , 



sup 

2(e 2 r)-"<t<T 



(t) ri \R A (t)\ < C ie -^ s sup \E(t) 



0<r<(e 2 r)- a 

+ CP 1 sup ((tP\E(t)\), 

(e 2 r)~ a <T<T 

sup (t) ri \R A {t)\ < D ( £ 2 r )-«(n+i) sup 

0<T<2(e 2 r)" a 



0<t<2(e 2 r)" Q 

Proof. To prove ( |5.5|) we begin with flPlD . Let 

A(t) =e-^ p{s)ds A(t). 

Then, A satisfies the equation 



d t A 



-e 2 TA + E{t) 



Solving (|5.8|) we get 



A(t) 



E(t) = e-ti p{s)ds E(t). 



e- £2rt i(0) + f e- £2T{t ~ s) E(s) ds 
Jo 



= e 



- £ 2 rt 



1(0) + R A (t). 



(5.5) 
(5.6) 

(5.7) 



(5.8) 
(5.9) 



(5.10) 
(5.11) 
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Below, in Proposition |5.2j we show that the real part of the integral of p{t) is uniformly 
bounded and of order 0(e 2 \ \ \W \ | | 2 ), for t > 0. Therefore, for some C > 0, we have by (|5.7|) and 



(5.12) 
(5.13) 



C- l \A(t)\ < \A(t)\ < C\A(t)\ 
C- l \E(t)\ < \E{t)\ < C\E(t)\ 

Consequently, it is sufficient to estimate A(t), in terms of E(t). 



Remark 5.1. Estimates of R a (t), which appears in the statement of Theorem |2.3| , are related 
to those for i?^(t) via : 

R a (t) = e - iXot+ ti p{s) ds R A (t) - (l - e R /o P^ ds ^j e - £2rt a(0). (5.14) 

Hence, by Proposition \5l\ , 

\Ra(t)\ < C \R A (t)\ + 0(e 2 \\\W\\\ 2 ) (5.15) 



From (|5.10| ) we have for any M > 0: 

rM rt 

lift) I < |A(0)|e- £2r *+ / e- £2rit - s) \E(s)\ds+ / e- £2r{t - s) \E{s)\ ds 

Jo Jm 

= \A(0)\e- £ ' 2rt + h(t) + h{t). 



Set 



M=(e 2 r)" Q , a>l. 



We now estimate the terms hit) and I 2 {t) in for 2{e 2 Y)- a <t<T. 

(t) ri h(t) = (t) ri [ M e- £2r{t - s) \E(s)\ds 

Jo 

< {t} ri e-^ £2rt ■ f M e -z 2 n\ts) ds . sup |^ (r) 

JO 0<T<(e 2 T)- a 



< sup ((t) ri e-^ 2rt ) • C^T)- 1 • sup \E{t) 

2{e 2 T)- a <t<T V 7 0<r<(e 2 r)- Q 

< Ce^ £2v) ~ S sup |£(r)|, 

0<r<(e 2 r)~ a 



(5.16) 



(5.17) 



for some 5 > 0. Therefore, 



sup «t) ri h(t)) < Ce~^ ' sup |£(r) 

2{e 2 T)- a <t<T 0<T<{e 2 T)- OL 



(5.18) 
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(t) ri h(t) < (t) ri f e - £2v{t - s) {s)- ri ds sup ((r) ri £(r)) (5.19) 



We estimate on the interval 2{e 2 T) a < t < T as follows: 

rt 

'(e 2 r)-<* ' ' (e 2 r)-«<r<T 

The integral is now bounded above using the estimate 

rt 



(t) ri ( e" e2r(i - s) {s)~ n ds < C(e 2 T)-\ t > 2(e 2 T)- a . (5.20) 

J(e 2 T)- a 

This gives 



sup (t) ri I 2 (t) < C(e 2 T)- 1 sup (r) ri £(r) (5.21) 

2(e 2 r)- Q <t<T (e 2 r)- Q <r<T 



Assembling the estimates ( |5.18| ) and ( |5.21| ) yields estimate ( |5.5| ) of Proposition |5.1| provided 



that ( [5 . 1 2| ) and ( p,13| ) hold. Estimate ( |5.6| ) is a simple consequence of the definition of -Ryi(t). 

Thus it remains to prove ( |5.12j ) and ( |5.13j ). By Q5.7| ) and Q5.9| ) it is necessary and sufficient 
to verify the following proposition: 



Proposition 5.2. Assume hypotheses (H1)-(H6). If p is given by fl4.12p then 



3? /' p(s) rfs < C^IHJyill 2 , t > 0, (5.22) 
Jo 

for some constant C depending on C, r± and £; see (H6). 



Proof of Proposition |5.2j: Using the estimates Q8.7|) and ( |8.9| ) we can infer that, p(t), given by 



Q4.12| ) is a series which converges uniformly on any compact subset of 1R. For each fixed t, it 



can therefore be integrated term by term to give 

rt r 2 rt 

E 



3? 



f p(s)ds = i fe^-^^MoiiHo-Xo-Pj-iOr'PcP^ds 



T E -Z" i (^o,(^o-A -/i,-^0)- 1 P c /3^ ) (5.23) 

4 j,fcez,jyfc M fc Mi 



Define 



pj,k = [Mo, {H - A - /ij- - iO) ^c/Jj^o) • 

/^fc — Pj 



Then ( |5.23|) can be expressed as: 

.2 ,.2 



$[ p (s)ds= £ - £ »ft.* = ir E ^fc + Pfc,i)- (5.24) 
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Now, since 



Pk,. 



-i(Hk—Hj)t _ I 
/If. — /Ij 



•■00, (H - A - p k - *0) -^PcPkiio 



Pk — A*j 



(/? fc Vo, (#o - A - /ifc + iO) i-PcPjipo, 



we have 



Pk ~ Pj 

Moreover, by flS.5| ) we can infer 

^(p ilfc + Pkj) = H (e*^-^)* - l) ^ fc + 23 ( e -*^-«)* - l) 5 j>fc , 
where, for j ^ k £ 2L, 

\ - p k - zO)- 1 P c /3 j o ) 
and for j ^ k, j £ 2Z, G 7 re s, 



fc^O, (-Ho - A - //j - ZO) 1 - (i/o - A - /Ufc + iO) ^c/^o 

(5.25) 



Pi,fc = (A^o, (-Ho - A - Pj - iO) 1 - (if - A c 



7T 



Pfc — Pj 



(Pji^o, S{H - A - Pk)Mo) • 



(5.26) 



(5.27) 



Thus, by d^j) and ( ggg ) 



/ p(s)ds = - E ^ (e^*""'* - l) Phk + - E 53 (e-ib*-** - l) <5 i)fc . (5.28) 

We now derive a uniform bound for 9ft J p(s) c?s. 

Estimating the modulus of the above sum, we have for any t: 

^ f t p(s)ds <y E M + f E M- ( 5 -29) 

By (H6), 

E M < < 1 1 \W\ II 2 . (5.30) 

We now bound the first term in ( |5.29| ). This requires an estimate of: 
1 



I Pi,* 



Pk — pj 



(3 k ip 0j (H - A - Pj - iO) 1 - (i7 - A - Pfc - iO) ^c/^o 
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for j ^ k E TL. We rely on the hypothesis (H3b) (singular local decay estimate 
which imply smoothness of the resolvent of Ho near accumulation points in cr con t(-f^o) of the set 
{A + Pjjjez- 

In order to treat both Ao + fij E cr cont (i/o) an d Ao + p>j <7 C ont{Ho) case simultaneously we 
regularize pj tk : 



T] 

Pi,k = 



1 



(/3fcVo» ( H o - A - fij - if])' 1 - (H - A - Hk - ^)~ 1 P c/ 9j'0o) ■ 

Clearly p j>k = lim^o p], k 

Now by the standard resolvent formula we have: 

Pj,k = {Mo, (Hq - A - Hk ~ iTlY l (H - A - pj - n7) _1 P c /^o) • 
Thus, using the singular local decay estimate (H3b), we get: 

(w + f3 k ipo, w_e~ tH ° s (H - A - pj - z?7)~ 1 P c iu_iu + ,3 i Vo 



(5.31) 



\PjM 



lim 

v\o Jo 



< lim 

n\o Jo 



-rjs 



ds 



< Ww+pkWWw+Pj 



\w 



_e- iHoS {H - A - Pj - iO) -1 P c u;_||ds 



< C||w + /3jfc|| / (s) ri ds 

J 

for some constant C depending on C and r%. Summing on j, k ElL, j ^ k yields: 

E \PiA < c\\\w\\\ 2 , 

j.fcez, j^k 

for some C > 0; see (|2.7|). Use of the bounds (|5.30|) and ( |5.33| ) in (|5.29|) gives 



(5.32) 



(5.33) 



K P pis) ds 
Jo 



< C e 2 \\\W\ 



for some constant C depending on C, r± and £. 



This completes the proof of Proposition |5.2j and therewith Proposition 5.1 
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5.2. Estimates for A(t) under the hypotheses of Theorem |2.1 



In this subsection we work under the hypotheses of Theorem |2.1| . In particular, we drop hy- 
pothesis (H6). We shall reuse the notation A and E for functions which are different from but 
related to those defined in section 5.1. 



Proposition 5.3. Suppose (H1)-(H5) hold. Then A{t), the solution of ( \4.11\ ), can be ex- 
panded as: 



A(t) = efo«M ds (e-^-fo^^AiO) + R A (t) 

R A (t) = f e-^t-rj-f* 3 />(»)«& e(t) dr, 
Jo 



where 



a(t) = ~U £ e*^-«)* ( Mo, S(H - A - ^)&Vo) ■ 



(5.34) 
(5.35) 

(5.36) 



is a real almost periodic function with mean M(o~) = 0, 7] in fl5.4ffi ) is a function whose real part 
has a bounded time-integral of order 0(e 2 \ \ \W\ | | 2 ) and E(t) is given in see also ( \4.1$) . 

For any a > 1, there exists 5 > such that R A (t) satisfies the estimates: 



sup (t) r i\R A (t)\ < C ie -( £2r/2) ~ 6 sup \E(t)\ 

2(e 2 V/2)~ a <t<T 0<r<(e 2 r/2)- a 

+ C(^r)- 1 sup ((r) n \E(r)\) 

(e 2 T/2)- a <r<T 



sup (t) ri \R A (t)\ < D (s 2 T/2)- a{ri+ V sup 

0<<<2(e 2 r/2)- a 0<r<2(e 2 r/2)- 

Proof. As in the previous subsection we begin with the equation for A(t): 

d t A(t) = (p(t) - e 2 r)A(t)+E(t), 



\E(t) 



(5.37) 
(5.38) 

(5.39) 



where p(t) and E(t) are given by ( [4.12| - f4"T3"D . In the previous section we transformed away the 
term p{t)A{t) using the "integration factor": exp(/Q p(s) ds). Under the current hypotheses, 
this can't be done because without (H6) 3? Jq p(s) ds may be unbounded as t — > oo, which could 
cause the estimates (|5.12| - |5~T3|) to break down. Instead, we proceed by a more refined analysis 
of p(t), which we now outline. 

We express p(t) as p(t) = e 2 a(t) + i](t), where n(t) has a time integral whose real part can 
be bounded by the estimates of section 5.1 and a part, e 2 a(t) which is almost periodic and of 
mean zero. Using this decomposition of p(t) we write ( |5.39| ) as 



d t A(t) 



-e 2 T + e 2 a{t) + V (t) A(t) + E(t) 
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Next introduce the change of variables 

A(t) = e-fi^ ds A(t) 



and obtain a reduction to 



d t A 



-e 2 T + e 2 a(t) A + E(t) 



E(t) = e-ti^ s)ds E(t). 



(5.40) 

(5.41) 
(5.42) 



With this strategy in mind we now proceed to derive the decomposition of p(t). We are 
mostly interested in its real part, so we start with it. 



4 



,ipo, (H - A - ft - iO) ^cPjipo 



Z e*^-«)* (foVo, (#o - A - m - z0)" 1 P c /3^ ( 



e 2 

= ~T E %Vj,k 
e 2 

= -o- Z ^(Vj,k + Vk,j)- 

° j,k&Z,jjtk 



In a manner similar to the derivation of ( |5.25| ) from ( 5.24 ) we find 

= Qe^-H)* (p h ^, (H - A - \i h + zO)- x P c /^o 
Using QB.5| ) in ( |5.44j ) and then replacing it in ( |5.43 ) we get 



-s is 



E 



(5.43) 



(5.44) 



j,kez,j^k 

_2. 



k i/j 0} \(H - Aq - /ij - iO) 1 - (# - A Q - fi k - iQ) x j P c (3j^ 



U 7](t) +e 2 a(t). 



Therefore, 



p(t) = ftp(t) + z9p(t) 
= ?7(t) + eV(t) 



(5.45) 
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where 



T](t) 



a(t) 



i 3 p{t) - -e 2 % e^-^ 1 



• k ij , [(H - A - Pj - 1O)- 1 - (H Q -Xo-Hk- iO^Pc] /3,'Vo) (5.46) 
% £ e i&*-«)* ( ^ 5 (# _ a - to)Pji>o) , 



see also (|5.36|) . 

Note that 9ft f r}(s)ds is uniformly bounded in t. To see this, recall the definition of pj^ in 
Lemma [572] (see ( |5.26| )): 



Pi.fe = (Mo, \(Ho -Xq- Pj -iO) 1 - (H - A - Pk - iO) x ] PcPjipo 

Pk — Pj v L J 

By (877), 3ft 77 (i) given by (|5.46|) , converges uniformly on t G IR. Therefore, for each t G 1R we 
may integrate the series term by term to obtain 



f ' rj(s)ds = -e 2 Y SR(e^-^*-l 

jo 8 ,rf,, v 



Pj,k- 



(5.47) 



Moreover the modulus of the right hand side in ( |5.47| ) is less or equal than \e 2 J2j.k j^k \Pj,k\ 
which by ( |5.32| ) is bounded by Ce 2 |||iy||| 2 for some constant C depending only on C and r\. 
Note that we derived ( |5.32j ) by using only hypothesis (H3b) and not relying on (H6). 
Thus we have 

& / r)(s)ds < Ce 2 |||iy||| 2 . (5.48) 
Jo 

To summarize, we have split p(t) into 

p(t)=r)(t) + e 2 a(t), 

such that ( |5.48|) is valid. If we now define A as in ( |5.40|) then, by (|4.11|) A satisfies the equation 
( 15.411) . Solving (|5.41| ) we get 

Jo 

: e - £2rt+£2 fo^ ds A(0) + R A (t). 



A(t) 



(5.49) 



From ( 5.42 ) and ( |5.4§ ) it is sufficient to estimate Ra(^)-, i n terms of E(t). 



Remark 5.2. The estimates of R a (t) which appears in the statement of Theorem |2..2| , are 
related to those for irU(i) via: 



Ra(t) 



- iX ° t+ Io r '^ ds R A (t) + (l - e ^(/ ^W<fa-7*)+«/ \W<iA e -e 2 (r-7)* a ( ). (5.50) 
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Before we estimate Rj\(t), we review some properties of the function cr(t). 

a(t) is an almost periodic function since the sum of the moduli of its Fourier coefficients is 
finite. Namely, by ( |2.20|) , the terms in the series ( |5.36|) defining a(t) are majorized by those 
of a convergent series (whose sum is CV -1 !!!^!!! 2 ). Therefore, the series in ( |5.36| ) is uniformly 
convergent. As the uniform limit of almost periodic functions, a(t) is then itself almost periodic, 
bounded by 

sup|a(t)| < C|||W||| 2 (5.51) 
ten 

for some constant C; see also section ^|. Moreover, a(t) has mean value zero since all the Fourier 
exponents are nonzero; see ( |5.36| ) and section |9|. Therefore 

r f r 

/ <j(s)ds<-(t-r), for t-r > M (5.52) 

provided M. is taken sufficiently large. It can be shown (see section |9] or @, page 42) that ( |5.52j ) 
holds provided 

M > "^'lf^». (5.53) 

r/2 

where L(T/4) (see Definition |9.1|) is such that in each interval of length L(T/4) there is at least 
one r/4-almost period for a . 

Using ( |5.51| ) and then (H5), we can choose: 

M = 8CL{T/A)/0 O (5.54) 

independently of e and still satisfy ( |5.53|) . 

We now return to the estimation of Ra- We split the integral in fl5.35| ) into two integrals, 
one from to t — M. and the other from t — A4 to t. For the former we use (|5.52j ) while for the 
latter we use ( |5.51| ). The result is 

\R A (t)\ < / e-^ v ^\E(r)\dr 
Jo 

+ f e £2 ( c IH w lll 2 - r )^|E(r)|rfr. (5.55) 

Jt-M 

The first integral in ( 5.55| ) can be bounded exactly as the term J * e~ £2v ^~ T ^\E{T)\dT in the proof 
of Proposition O. The second integral in ( |5.55| ) is bounded in the following manner: 

u\n f e - 2 ^lll^lll 2 -r)(*-)|E(r)|dr< 

Jt-M 

77%fcf e^ w ^^dr sup «r>^|^(r)|) 
{t — Jviy 1 Jt-M t-M<r<t 

< D sup ((r) ri |£(r)|). (5.56) 

(e 2 r/2)- Q <r<t 
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Note that e and consequently e 2 T ~ £ 2 |||VT||| 2 are small, so we can consider J\4 <C (e 2 T /2) Q 
and D ^ M. <^ (e 2 r) -1 . The result is ( |5.37|) . A simple bound, using the definition of iiU(t) 
yields ( pg ). 

This completes the proof of Proposition |5.3j. 



6. Dispersive Estimates and Local Decay. 

In this section we prove the local decay of (f>d and the decay in time of the remainder terms, 



E(t), in bound state amplitude equation ( 4.11 ) of section |j. The arguments rely on hypotheses 
(H1)-(H5) and results of the previous section, so we will handle Theorem [Tl] first. However, 
due to the differences between Theorems [2.2| and |2.3| we separately finish their proofs in the final 
two subsections of this section. We will repeatedly use the following: 



Lemma 6.1. For any rj G [0, r\] and j G Zj we have 

w_e- lHo{t - s) P c f{s)ds 

Jo 



< C(t)- r > sup «t)"|K/(t) 

0<T<t 



(6.1) 



and 



f w_e- iHo{t ~ s) P c {H - A - fij - «0)- 1 /(s) ds 

J 



< C(t)-*> sup «T>'||u; + /(r)||). (6.2) 

0<T<t 



Proof. The proof follows from the assumed local decay estimates on e lHot ; see (H3a). Namely, 
using that r± > 1, 



f w ^e~ lHoit - s) PJ(s) ds 
Jo 



< 


L 







• sup «t)»»|K/(t)||) 

0<r<t 

< C [\t-s)-^(s)-v ds sup ((rr\\w + f(r) ||) 

JO 0<T<t 

< C(t)-* sup ((rr\\w + f(r)\\) 

0<T<t 

which proves (|6.1|) . The proof of ( |6.2|) is identical, and uses the singular local decay estimate of 
(H3) (b) []. 

We now define the norms 

(6.3) 



[A} a (T) = sup (r} a \A(r) 

0<t<T 



and 



Then we have 



l >d\LD, a ( T ) = SU P ^TWw-^dir) 

0<t<T 



(6.4) 
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Proposition 6.1. For any T > and n G [0, r±], 

IMld, v ( t ) <C( \\w+(f) d (0)\\ + \ £ \ \\\ W \W i A U T ))- ( 6 - 5 ) 
Proof. From equation (|4.9|) we get, using the assumed local decay estimate for e~ lHot and (|6.1|) , 



\w-Mt)\\ < E IK<M*)!I 

j=0 

< C(t)-*>\\w+(j> d (0)\\+ C\e\(t)-* [A] v (t) sup \\w + W(s)^ \\ 



0<s<t 



+ c Hill w\\\ (t)-*[MLD, n ®- 



(6.6) 



Since ||u> + W(s)'0o|| < 1 1 1^1 1 1 HV'oll = lll^lll an d \s\ |||W^||| is assumed to be small, mul- 
tiplying both sides of this last equation by (t)' n and taking supremum over t < T yields ( |6.5|) . 



We now estimate E(t). 
Proposition 6.2. Let T > 0. For any 77 G [0, ri\: 

[E] V (T) < C(e 2 \\\W\\\ 2 \A(0)\ + \s\ \\\W\\\ |K0 d (O)|| + \e\ 3 \\\W\\\ 3 [A^T) ) (6.7) 
Proof. 



E(t) is defined in ( [4. 13| ) . From these equations it is seen that we need to bound the following 
terms: 



Ri = \e 2 \A(0)\ E \(Pk4>o, e~ mot (H - X - fij - iO^P^o) 

4 j.A-ez 



and 



|e(^o, H/(t)0 2 )| 



e(W(OVo,e-^°V rf (0)) 
e 2 (W(t)ip , f e- iHo ^P c W(s)Ms)ds 



The estimates of the above terms repeatedly use Lemma |6~T|. Let 77 G [0,ri]. 
Estimation of Rj: 
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R x = -e 2 \A^)\ E I ( w +Mo, w_e- lHot (H - A - m - zO)- 1 P c u;_ w+Ptfo 



2\\\rxr\ 1 1 2 



< C7|A(0)|^|||W|||^(tp 



(6. 



by the local decay estimates 
Estimation of R2 



From ( 4.11 ) we have that 



\d s A(s)\<C\e\ 



\A(s)\ + \E(s) 



(6.9) 



since is linear in \e\ 



and 3f?p, T are quadratic. 



Applying Lemma 6.1 to ^2 we then get 



7^ E 



w + P k ipo, / w_e 



-iH (t-s) 



(Hq — \q — — iO) 1 P c W-d s A(s)w + Pjijjods 



<*>-"(kl lll^lll M„(*) + [#],(*))■ 



(6.10) 



Estimation of |e (^ , W (t)^ (t))\: 



Since, by definition, 4>a(0) = P c 0d(O) we can apply local decay estimates for e %Hot to get 
|e(W(t)^, 0o(*))| < C\e\ \\\W\\\ (t)^ |K0 rf (O)||. (6.11) 



Estimation of |e (V>o, W(t)02)| 



(6.12) 



Applying Lemma |6J] as before we get, for < t < T, 

W(t)<h)\ <Ce 2 \\\W\\\ 2 (tV [Mld, v (T). 
Using Proposition |6.1| to estimate [4>d\LD,r)(t) in (|6.12|) , we get 

| £ ^o, W(t)<h)\ <Ce 2 \\\W\\\ 2 ||«> + ^(0)|| + |e| |||W||| [A] v (t) } . (6.13) 

Finally, combining the above estimates, we can bound [E] V (T) for any 77 G [0, ri] as follows: 



[E] V (T) < C{e 2 \\\W\\\ 2 \A(0)\ + \e\ \\\W\\\ |Kfc(0)|| + e 2 |||iy||| 2 [E]„(T) + |e| 
Since |e| |||W||| is assumed to be small, Proposition |6.2| follows. [| 



3 [A] V (T)} 
(6.14) 
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We can now complete the proof of Theorem |2.1| . To prove the assertions concerning the 
infinite time behavior, the key is to establish local decay of <pd, in particular, the uniform 
boundedness of [4>d\LD,niT). This will follow directly from Proposition |D] if we prove the 
uniform boundedness [Aj ri (T), or equivalently [A] n (T). 



Proposition 6.3. Under the hypothesis of Theorem 2.1, there exists an Eq > such that for 



each real number e, \e\ < Eq there is a constant C#, with the property that for any T > 

[A] n (T)<C* 

Proof. 



We begin with the expansion of A(t) given in Proposition |5.3| . Multiplying ( 5.34Q by (t) 
and taking the supremum over < t < T we have: 

[A] n (T) < c(\A(0)\ (e 2 T/2)^ + sup ( t )^\R a (t)\ + sup (t)*\R A (t)\ 

V 0<r<2(e 2 r/2)- a 2(e 2 r/2)- Q <r<T 



ri 



(6.15) 



The right hand side of (|6.15|) is estimated using Proposition |5~3 . 



[A] ri (T) < C\A{0)\ {E 2 T/2)- ri +D (E 2 T/2)- a ^-^ [E] (2(E 2 T/2)- a ) 
+ d e-^M- s \EU2{E 2 T/2)- a ) +C 2 (E 2 T/2)- 1 [E] ri (T). 

Next, we apply Proposition |6.2| which yields: 

[A] ri (T) < C\A(0)\ (e 2 r/2)- ri + D ( £ 2 r/2)- Q(ri+1) [E] Q (2(e 2 T /2)~ a ) 

+ C ie - (£2r/2) ^ [E} (2(E 2 T/2y a ) (6.16) 

+ ^(^r^)- 1 ^ 2 ^^)!!!!^!!! 2 + | e | |||w||| |K<Mo)|| + k| 3 |||iy||| 3 [A] ri (T)). 



Note that by Proposition [T2] and the simple bound: 

[A]o(T) < Hoi 

[E]o(2(e 2 T /2)~ a ) is bounded in terms of the initial data and |e| \\\W\ 
Choose £ such that: 

i "p - s = o, 

where C 2 is the same as in ( |6.16| ). Then, for |e| < eq 



[A] ri (T)<a (6.17) 
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Here, C* depends on \\4>o\\, \\w + (f>o\\, r\ and e. 



This completes the proof of Proposition |6.3| and therewith the t — > oo asymptotics asserted 
in Theorems |2.1| - |2.3| . 1 1 



It remains to finish the proofs of the Theorems and |2.2| . Due to some differences we 
consider them separately in the following two subsections. 



6.1. Proof of Theorem 2.3 



In order to obtain ( |2.23| ) we note that ( |4.7p , ( |5.3| ) and ( |5.14| ) together with the definition of oo(t) 
in ( |2.16| ) already gives us: 



a(t) = e - iAof +/o P^ ds (A(0)t 



-e 2 Ft 



+ R A {t 



a(0)e- £2rt e iuj{t) + R a (t). 



which is in fact the second relation in (|2.23|) . The third is a direct consequence of the second 
since P(t) = |a(t)| 2 while the fourth relation is exactly ( |4.9| ). 

It remains to prove the intermediate time estimate (|2.17|) . The ingredients are contained in 
and its proof. First, by (|5.15| ) 



\R a (t)\ <C\R A {t)\ + 0(e 2 \\\W\\\ 2 ) 

So, it suffices to prove an 0(\e\ |||W|||) upper bound for Ra- 
Using (|5.4j ) we know that 



\R A (t)\ < I e- e2T{t -^ \E(t)\ dr. 
Jo 



(6.18) 



Let T Q denote an arbitrary fixed positive number. We estimate (|6.18 ) for t G [0, T (e 2 r) _1 ]. We 
bound the exponential in the integrand by one (explicit integration would give something of 
order (e 2 r) -1 ), and bound |-E(t)| by estimating the expressions in the proof of Proposition |6.2| . 
First, the estimates of Proposition ^72] for Ri and \e (ipo, W{t)cf) {t)) \ are useful as is. Integration 
of the bounds ( |6.8| ) and ( |6.11j ) gives: 

/* e- £2r{t ~ T) R 1 dr < C e 2 \\\W\\\ 2 \\w+t 
Jo 



-e 2 T(t-r) 



e(ip , W{t)(f> {t))\ dr < C \e\ \\\W\\\ \\w+ 



(6.19) 



To estimate the contributions of i?2, first observe that by ( |6.9|) and Proposition |6.2| with r\ = 

\d s A(s)\ < C \e\ \\\W\\\ ||w+0(O)|| (6.20) 
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Therefore, using local decay estimates we have: 

C e- e2r ^R 2 dr < C To^r)" 1 \e\ 3 \\\W\ 
Jo 



< D\e\ \\\W\\\ \\w. 



Finally, we come to the contribution of \e (ipo, W{t)(p 2 )\. We rewrite it as follows. 



\e(ipo, W{t)cj> 2 )\ 



(W(s)e iHo ^P c W(t)ijj , d (s)) ds 

f ' e 2 (w + W(s)w + ■ w_e iHo{t - s) P c w. ■ w + W{t)^ , w_<p d {s)) dsk.21) 
Jo 



Recall that by ( |4.9|) (p d = (f) + <pi + <p 2 , where (fio(t) = e~* H °*0d(O). Using local decay estimates 
(H3a), the contribution of the term (fio(t) can be bounded by C e 2 \\\W\\\ 2 ||u> + 0d(O) || (r)~ n . 
Multiplication of this bound by e~ £2r ^~ T ^ and integration with respect to t gives the bound 
C e 2 |||W / ||| 2 ||w + ( i(O)||. To assess the contributions from (pi + 4> 2 , note that local decay estimates 
(H3a) imply 

||w_(0i + 2 )|| < C \e\ \\\W\\\ \\w+<f>(0)\\. (6.22) 
Putting together the contributions from <ft and from (f>i + <f) 2 , we have: 

e- £2vit - T) \e^ W(t)fa)\ dr < C ( e 2 \ \ \W\ \ \ 2 \\w+<l> d (0)\\ + (e'vy 1 \e\ 3 \\\W\\\ 3 ) 

(6.23) 

The above estimates and ( |5.15| ) imply ( J2.17 ) . Now, ( |2.1^ ) is a direct consequence of ( |2.17p and 
the relation P(t) = |a(t)| 2 . This concludes the proof of Theorem 



6.2. Proof of Theorem I2T2" 



As in the proof of Theorem |2.3| relations ( |4.7|) , ( 5.34Q , Q5.50| ) and the definition of u>(t) in Q2.16| ) 
gives: 

a(0)e" £2 ( r -^e^ w + R a (t), 



a(t) 



which is the second relation in ( |2.15| ). In what follows, the only difference from the previous 
argument is in estimating R a (t). 
We start with the relation (|5.50|) : 



Ra(t) 



- iXot+ Io^ ds R A (t) + fl - (fHJo^ds-^+StJ^v^ds^ e - £ *(r- 7 )t a (Qy (g 24 ) 
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Since a(t) is an almost periodic function with zero mean, for any 7 > there is an M.^ > 
such that whenever \t\ > A^ 7 

/ a(s)ds < jt. 
Jo 

On the other hand for \t\ < J\4 y , using (|5.51|) we have: 

[ t a(s)ds < C.M 7 |||W||| 2 . 
Jo 

So, in both cases, 

/ a(a)ds-jt < CMJ\\W\\\ 2 . 
Jo 

Substituting now in ( |6.24|) and tacking into account that by (|5.48 



3? (\(s)ds < Ce 2 \\\W\ 
Jo 

uniformly in t, we get 



Ra(t)\ < C\R A (t)\ + 0(e 2 \\\W\\\ 2 ) (6.25) 



It remains to prove an 0(\e\ |||W|||) for Ra^)- Looking now at (|5.55 ) we see that we can 



bound the exponential by max{l, e 52 ^"' 1 ^'" 2-1 "^}. Now, the same argument as in the end of 
the previous subsection will give us the required result. 
This completes the proof of Theorem [2.2| . 

7. Generalizations 

In the previous sections we considered perturbations of the form eW(t), with W(t) independent 
of e. In this section, we shall extend our theory to a more general class of potentials, W e , which 
are small for small e, but which may deform nontrivially as e varies. 
Consider a family of perturbations W and the general system 

id t <f>(t) = (H + W(t)))<Kt), 
<P\t=o = 0(0). (7.1) 

where W G W (compare to (|2.1|) . The results are 

Theorem 7.1. Suppose that H and any W G W satisfy hypotheses (H1)-(H5). In addition 
assume: 
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(H7) Equi-almost periodicity: There exists a positive constant Lg Q , independent of W G W, 
such that in any interval of real numbers of length Lq , the function | \\W\ | |~ 2 cr(t) (\ \\W\\ \ ^ 0), 
where 

a{t) = -j^ e^-^iMo^Ho-Xo-fji^o). (7.2) 

has a 6q/ A almost period, 60 is given by (H5). More precisely, there exists Lg > which does 
not depend on W such that in any interval of length Lq there is a number r = r(# /4), such 
that for all t G IR 

\\\W\\\- 2 a(t + T) - |||W|||" 2 <7(t) I < o /4. (7.3) 
If w + (p(0) G H, then, there exists an Eq > (depending on C, r±, 9 and Le ) such that 



whenever |||W||| < Sq, the solution of ( |7. 1\ ) satisfies the local decay estimate (identical with the 
one in Theorem 12.11 ): 

Ik'- (f>{t)\\ < C(t}- ri \\w+ 0o ||, t G IR. (7.4) 

Sketch of Proof. Once we drop e from all expressions (since it is not present in the actual setting), 
the arguments in the previous sections hold in this case except the analysis of a(t) in Proposition 
1| Formulas ( |5.36| ) and ( |7.2j ) are the same, but now fj,j, j3j, j G TL are not fixed as they define 



W by (H4) and W sweeps a general class W. This may prevent us to find a fixed time interval, 
M., independent of W G W, after which a(t) is within T/2 distance from its mean, see relations 

Nevertheless, (H7) is exactly what we need to overcome the difficulty. A straightforward 
calculation shows that any 8q/4 almost period of 1 1 \W \ \ \~ 2 a(t) is a T/4 almost period for u{t). 
Consequently, L(r/4) in (|5.54[) is bounded above by L(6q) given in (H7). But the latter is 



fixed, so, we can choose 

M = 8CL(9 )/9 ( 7.5) 

independent of W G W and still satisfy ( |5.53| ) hence ( [5.521 ). 

Finally, we can close the arguments exactly as we did for Theorem |2.1| . 



Remark 7.1. Theorems analogous to Theorem |2.2| (respectively Theorem |2.3| ) can be proved 
under hypotheses (H1)-(H5), (H7) (respectively (H1)-(H6)| 

Examples: (H7) holds trivially for 
(1) W = {eW(t,x) : e G IR, W fixed} 
or 



39 



(2) W = {eW(£-H, x): e e JR - {0}, |e| < 1, W fixed}. 

In the Example 1, e cancels in the formula 1 1 \W\ \ \~ 2 a(t) while in the Example 2 we have a 
time dilatation which shrinks the gaps between the almost periods, so L(9 ) valid for W is good 
for the entire family. 

(3) There are more general families of perturbations W for which (H7) holds. For example if 
W is equi-almost periodic, see section []. 

8. Appendix: Singular operators 

In this section we present the definition and the properties we needed previously for the singular 
operators: 

e -iH t ( Hq _ A _ . Q) -l 6 ( Hq _ A) p V ( Hq _ A) -l 

and establish the identities 

(Ho-ATiOy^c = P.V.(H -A)- 1 P c ±m5(H -A)P c 
suggested by the well known distributional identities 

(x^iO)- 1 = P.V. - ±m 8(x). 

x 

Recall that we are in the complex Hilbert space 7i with self-adjoint "weights" w± and 
projection operator P c satisfying (i), (ii) and (iii). We can then construct the complex Hilbert 
space H + as the closure of the domain of w + under the scalar product (f,g) + = (w + f,w + g) 
and the complex Hilbert space 7i_ as the closure of P C TC under the scalar product (f,g)_ = 
(w-f,w_g) . 

By hypotheses of section ||, H is a self-adjoint operator on 7i and satisfies the local decay 
estimate fl2.3|) . Based on this property, in jLTJ ^2], ^3f it is proved that for A in the continuous 
spectrum of H and t G IR 

T t = i\im e- i{H °- A - iv)s dsP c 
v\o Jt 

TT = -i lim /"* e- i{H °- A+iri)s dsP c 

are well defined, linear bounded operator from 7i + to We then define 

e- iHot (H - A-i0y l P c = e- iAt T t (8.1) 
e +^ot( jH - o _ A + i0 )-ip c = e +iAt T t *, (8.2) 
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and 

P.V.(iJo-A)- 1 P c = ^(T + T*) (8.3) 

5(H -A)P C ee ^-(T -T*) (8.4) 

Note that the definitions imply the identities: 

(Hq - A T zOr 1 P c = P.V. (Hq - A)' 1 P c ± tn5 (H - A) P c . (8.5) 

Particularly important properties of these operators are their symmetries when viewed as quadratic 
forms on 7i + x H + . For example, on any f,gE H + the quadratic form induced by T t is given 
by: 

(/>#) (w+f,W-T t g). 

Note that 

lim (/, T?g) ee lim f /, z /°° e-^- A -"^ S P c s) = (u;+/, «;_T t (/) (8.6) 

»?\0 r?\0 \ Jt J 

by the following calculation: 

lim (f,i r e- i{H °- A - iri)s dsP c g) = lim ( f, P c i f°° e - i{Ho ~ A ~ iri)s dsP c g) 
v\o \ Jt J ti\o V Jt J 



lim (f,w+w-P c i [°° e~ i{H °- A - iv)s dsP c g) 
v\o \ Jt J 

lim (w + f,w„i r e- i{Ho - A - ir ' )s dsP c g 
v\o \ Jt 



= (w+f,w-T t g), 

where we used the that P c is a projection operator commuting with the integral operator, the 
identity w + W-P c = P c on TC, the self adjointness of w± and P c , and lim^\o W-Tl' = W-T t in 
£(H + ,H). 

Identity (|8.6| ) suggests the notation: 

(f,g)^(f,T t g) 

for the quadratic form induced by T t , where (•, •) can formally be treated as the scalar product 
in 7i. Moreover, ( |S.6|) implies 

(f,T t g) = (T;f,g). 

Therefore, the quadratic form induced by P.V.(H — A) _1 P C is the symmetric part of the one 
induced by T while 5(H — A)P C induces the skew-symmetric part of it divided by the factor 
in. As a consequence both the forms corresponding to the last two operators are symmetric. 
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In conclusion, for any /, g G Domain(w + ), t G IR and A G cr cont (if ) we have: 

< C t \\w + f\\ \\w + g\\ (8.7) 
(f,5(H -A)P c g) = (w + f,W-5(H -A)P c g) < ^ \\w + f\\ \\w + g\\ (8.8) 

7T 

(/,P.V.(fZ- -A) _1 Pc5) = (/^-P-V.^o-A)-^) < Co IK/H ||«; + ff||.(8.9) 

The inequalities are due to the boundedness of T t , where C t denotes the norm of T t in C{7i + , 7i ). 
Moreover, the following symmetry properties hold: 

(f^^iHo-A^iO)- 1 ?^) = (e ± * Hot (H -A±tO)- 1 PJ,g) 
(f,5(H -A)P c g) = (5(H -A)PJ,g) 
(/^.V.^o-A)-^) = (P.V.^o-Ar'Pc/,*?). 

9. Appendix: Almost periodic functions 

In this section we present the definition and the properties of almost periodic functions we used 
throughout this paper. We will confine to functions of the form / : IR — > X where X is a 
complex Banach space with norm denoted by || • ||. 

Definition 9.1. We say that 

f:JR^X 

is almost periodic if and only if it is continuous and for each e > there exists a length 
L(e, f) > such that in any closed interval of length greater or equal than L(e, f) there is at 
least one r with the property that for all t G IR we have 

\\f(t + r)-f(t)\\<e. (9.1) 

The number r with the property above is called an e-almost period for f . 

We say that the family, T of almost periodic functions is equi-almost periodic if L(e, f) can 
be choosen independently of f G T . 

Example 1 Any continuous periodic function is almost periodic since for any e > we can 
choose the length L(e) to be the period of the function. 

Theorem 9.1. Any almost periodic function has a relative compact image. 
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The proof of the theorem can be found in |], Property 1, pp. 2]. In particular, any almost 
periodic function / : IR — > X is in the Banach space of all bounded and continuous functions 
on IR with values in X, C(X), endowed with the uniform norm. The next result is Bochner's 
characterization of almost periodic functions, see for example ||, Bochner's theorem, pp. 4]. 

Theorem 9.2. (Bochner) Let f : IR — > X be a continuous function. For f to be almost periodic 
it is necessary and sufficient that the family of functions {f(t + h)}, — oo < h < oo is relatively 
compact in C(X). 

As a consequence of Bochner's criterion and Property 4 from pp. 3] we have: 

Theorem 9.3. Suppose X\, X2, . . . , Xk+i are Banach spaces, : IR — > Xi, 1 < i < k are 
almost periodic functions and g : Yii =1 — > Xk+i is continuous. Than g(fi(t), f2(t), ■ ■ ■ , fk(t)) is 
an almost periodic function. 

The last theorem has very important consequences in the theory of almost periodic functions. 
We will list only those which are useful in our presentation. 

Corollary 9.1. A finite sum of almost periodic functions with values in the same Banach space 
is an almost periodic function. 

Corollary 9.2. A product between a complex valued almost periodic function and an arbitrary 
almost periodic function is an almost periodic function. 

Corollary 9.3. If 7i is a complex Hilbert space, C(TC) is the Banach space of the bounded 
linear operators on TC and W : IR — » C(7i) is an almost periodic function then for any <p, if) G 7i 
the following functions are almost periodic: 

t -> W(t)tp 

t - w{t)<p) 

t ^ (W(t)iJ>, W{t)<p) , 
where (•, •) denotes the scalar product on 7i. 

Another essential result in the theory of almost periodic functions is (see for example || 
Property 3, pp.3]): 
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Theorem 9.4. Any uniform convergent sequence of almost periodic functions converges to- 
wards an almost periodic function. 



Corollary 9.4. If {/^}j G z C IR and {/^Jjez Q X satisfies S je z WfijW < °°> then 

E ^ 

is an X— valued almost periodic function of t. 

Proof: According to Weierstrass's Criterion the series X^ G z e l/ ^*/3j is uniformly convergent on 
IR. 

By Corollary |9.1| and Example 1 the partial sums of the above series are almost periodic. 
The result follows now from Theorem |9.4| . 

We continue with the harmonic analysis results for almost periodic functions. 

Theorem 9.5. (Mean Value) If f : IR — » X is almost periodic then the following limit exists 
and it is approached uniformly with respect to a G IR: 



Moreover, whenever 



we have 



for all a G IR. 



lim - f f(s)ds = M(f) G X. 
;-»-°o t J a 

4su Pse]R ||/( S )||L( £ /2,/) 



t > 



\\M{f)- l -f a+t f{s)ds\\<e, 

l Ja 



The proof of the mean Value Theorem in this form can be found in % pp. 39-44]. Note that 
although Bohr's book consider only complex valued almost periodic functions the proof can be 
carried on to Banach space valued functions by simply replacing the modulus by the norm and 
the Lebesque's integral for complex valued functions by the Bochner's integral. 

The results of the next theorem are presented in || Chapter 2]. 



Theorem 9.6. (Fundamental theorem) If f, g : IR — > X are almost periodic then: 
(a) for any ft G IR, 

hm^ f f(s)e-^ s ds = a(u,J) 
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exists and is non-zero for at most a denumerable set of fi's; if a(/x, /) ^ then a(/x, /) is called 
a Fourier coefficient for f while u, is called a Fourier exponent; 

(b) a(/i, /) = a(/i, g) for all fi G IR if and only if f = g; 

(c) let A(/) = {fi : a(fi,f ) 7^ 0} denote the set of Fourier exponents for f; then there is an 
ordering on A(/), A(/) = {//i, fi 2 , . . .} independent of the Fourier coefficients, such that for 
any e > there exist the numbers N(e) G IN, < k nj£ < 1, n G IN with the property that the 
trigonometric polynomial 

N(e) 

p e (t)= ^KM^f)^ 

n=l 

satisfies 

\\f(t) - P e (t)\\ <e for all ten. 
Moreover k n ^ £ can be choosen such that for any fixed n, lim e \ k nt£ = 1 . 

In this paper we use a less general result than the above Fundamental Theorem, namely: 

Corollary 9.5. If f(t) = E jG z where {^} je z C IR and E je z \\Pj\\ < oo, then A(/) = 

{/(jljez- o(A*j, /) = /3j, j G !Z, in particular if fij ^ 0, j E then M(f) = 0. Moreover, we can 
arbitrarily order A(/) and still have that for any e > there exists a natural number N(e) such 
that: 

j=N(e) 

||/(t)- £ e*"/y<e, 

j=-N(e) 

or, in other words, in this particular case the conclusion in part (c) of the Fundamental Theorem 
is valid even if we have an arbitrary order on A(/) and we choose kj jE = 1. 

Proof: By the Weierstrass criterion the series 

f(t)e~ ilJt = ^jto-^Pj 

is uniformly convergent on IR. So, when we compute a(/i, f) we can integrate term by term and 
therefore use the identities: 

lim - /* e~ lXs ds = 
t Jo 

to get the first part of the Corollary. The last part is a direct consequence of the fact that / is 
an absolute and uniform convergent series. [] 



if A^O 

1 if A = 
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